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Abstract 

The gravitational radiation originating from a compact binary system in cir- 
cular orbit is usually expressed as an infinite sum over radiative multipole 
moments. In a slow-motion approximation, each multipole moment is then 
expressed as a post-Newtonian expansion in powers of v/c, the ratio of the 
orbital velocity to the speed of light. The "bare multipole truncation" of 
the radiation consists in keeping only the leading-order (Newtonian) term in 
the post-Newtonian expansion of each moment, but summing over all the 
multipole moments. In the case of binary systems with small mass ratios, 
the bare multipole series was shown in a previous paper [Simone et al. 1997, 
Class. Quantum Grav. 14, 237] to converge for all values v/c < 2/e, where e 
is the base of natural logarithms. (These include all physically relevant values 
for circular inspiral.) In this paper, we extend the analysis to a "dressed mul- 
tipole truncation" of the radiation, in which the leading-order moments are 
corrected with terms of relative order (v/c) 2 (first post-Newtonian, or 1PN, 
terms) and (v/c) 3 (1.5PN terms). We find that the dressed multipole series 
converges also for all values v/c < 2/e, and that it coincides (within 1 %) with 
the numerically "exact" results for v/c < 0.2. Although the dressed multi- 
pole series converges, it is an unphysical approximation, and the issue of the 
convergence of the true post-Newtonian series remains uncertain. However, 
our analysis shows that an eventual failure of the true post-Newtonian series 
to converge cannot originate from summing over the Newtonian, 1PN, and 
1.5PN part of all the multipole moments. 
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I. INTRODUCTION 



Measurement of gravitational waves from inspiraling compact binaries (composed of neu- 
tron stars and/or black holes) by forthcoming laser interferometric detectors (such as LIGO 
|H] and VIRGO @) will rely heavily on theoretical templates obtained from approxima- 
tions to the exact laws of general relativity ||£|] . One possible method to calculate these 
templates, which has been pushed to a high degree of accuracy is the post-Newtonian 

approximation, which is based upon the assumption that the orbital motion is slow. Another 
approach is black-hole perturbation theory |||| , which is accurate only in the unrealistic sit- 
uation of systems with small mass ratios. Currently, the post-Newtonian approach appears 
to be the most promising for generating the templates. Unfortunately, however, there is 
considerable evidence that the post-Newtonian approximation has poor convergence fTOJIT 



if it converges at all. It is important to understand the reasons for the poor convergence. 
Although this issue is still an open one, some preliminary steps were taken in this direction 



in a previous paper [p~2 1 . 



The amount of energy radiated per unit time (denoted E) by a binary system in circular 
orbit can be expressed as an infinite sum over the multipole moments of the gravitational- 
wave field, which are related to the behaviour of the source. In a slow-motion approximation, 
each multipole moment is then expressed as a post-Newtonian expansion in powers of v/c, 
the ratio of the orbital velocity to the speed of light. Since higher-order moments come with 
a higher power of the orbital velocity, the post-Newtonian approximation to E consists of 
truncating the multipole sum at a given multipole order, and in calculating the contributing 
moments to the appropriate post-Newtonian order. We will refer to this expansion of E 
in powers of v/c as the "true post-Newtonian series". At present, only a few terms of this 
post-Newtonian series have been calculated and the issue of its converge is not at 

all clear. However, since post- Newtonian corrections can be computed for all the multipole 
moments, the simpler issue of the convergence of the multipole sum can be addressed. 

In reference [12[], the convergence of the multipole series was studied in the restricted 
context of a "bare multipole truncation", in which the multipole moments were calculated 
only to leading-order in a post-Newtonian expansion. We refer to these bare moments as 
"Newtonian", and the bare multipole truncation therefore consists of discarding all post- 
Newtonian corrections to each of the multipole moments. When an approximation to E is 
constructed by summing over all these bare moments, it is found that the series converges 
for all orbital velocities less than a critical value. This value depends on the mass ratio of the 
system; in the case of very small mass ratios, to which we specialize in this paper, the con- 
dition for convergence is t>/c < 2/e ~ 0.7358, where e is the base of natural logarithms. The 
convergent interval includes all physically relevant values for the orbital velocity, which must 
be smaller than its value at the innermost stable circular orbit, v- lsco /c = l/y/6 — 0.4082. 
Although the bare multipole truncation is an unphysical approximation which compares 
poorly with the numerically "exact" results [IIJ, this analysis reveals that an eventual fail- 



ure of the true post-Newtonian series to converge cannot originate from summing over the 
Newtonian part of all the multipole moments. 

In this paper, we extend the results of reference to a "dressed multipole truncation" , 
in which the leading-order moments are corrected with terms of relative order (v/c) 2 (first 
post-Newtonian, or 1PN, terms) and (v/c) 3 (1.5PN terms). This improved analysis relies 
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on a previous paper of ours , in which the equations of black-hole perturbation theory are 
integrated for an arbitrary perturbing stress-energy tensor T a/3 (the unperturbed spacetime 
is assumed to be Schwarzschild) . Applying the results of this paper to the particular case of 
a particle moving on a circular orbit returns the desired post-Newtonian corrections to the 



bare multipole moments derived in reference [13 



The dressed multipole truncation gives three separate series for E. The first is the bare 
multipole series, in which we sum over the Newtonian part of all the multipole moments. The 
second and third series incorporate the 1PN and 1.5PN corrections to the bare moments, 
respectively. We find that each of the series converges for v/c < 2/e, and we conclude 
that inclusion of these post-Newtonian corrections does not affect the convergence of the 
multipole sum. This shows that an eventual failure of the true post-Newtonian series to 
converge cannot originate from summing over the Newtonian, 1PN, and 1.5PN part of all 
the multipole moments. 

The following sections of the paper provide the details of the calculation. The dressed 
multipole truncation is introduced in Sec. |J, an d the convergence of the resulting series is 
established in Sec. [TTI|. Finally, concluding remarks are provided in Sec. [TV]. 



II. DRESSED MULTIPOLE TRUNCATION FOR GRAVITATIONAL RADIATION 



In this section, we extend the bare multipole calculation of reference [|T2| by calculating 
the 1PN and 1.5PN corrections to the multipole moments of the gravitational- wave field. 
Such corrections have been calculated for general sources, first by Blanchet |14| using post- 



Newtonian theory, and then by us || using black-hole perturbation theory. To carry out 
this calculation, we apply the results of reference || (hereafter referred to as LP) to the 
specific case of a point particle moving on a circular orbit around a Schwarzschild black 
hole. We denote by fi the mass of the particle, and by M the mass of the black hole. It is 
assumed that p/M <C 1. We also denote the orbital radius by r , and the angular velocity is 
n = (M/r 3 ) 1//2 . Finally, we define the orbital velocity by v = fir = (M/tq) 1 ^ 2 . Throughout 
we use units such that G = c = 1. 

The gravitational-wave luminosity is given by (TT 
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where Xi m (u) and Si m (u) are the mass and current multipole moments, respectively, while 
u denotes retarded time. (Because the orbit is fixed, E does not actually vary with u.) To 
evaluate E, we shall use equations (5.37)-(5.39) of LP, which give the Fourier transform of 
the moments, Xi m (uj) and Si m (uj). These equations incorporate post-Newtonian corrections 
of up to order v 3 , relative to the leading-order, Newtonian expressions. 

To evaluate the multipole moments, we must first compute the source terms, which 
are constructed from the stress-energy tensor of the orbiting particle ||16|| . Using equations 



(5.24)-(5.27) of LP, we find that the only nonzero source terms are p = RA, _ij = _i JA, 
ij = i J A, and t = qTA, where 
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and i J = — (-iJ). 

The second step is to substitute these results into equations (5.37)-(5.39) of LP. After a 
few steps of algebra, the mass and current multipole moments become 
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are the spherical harmonics of spin weight 1 and —1, re- 



while iYi m (9,(f) / 
spectively |TT|| . 

The third step is to invert the Fourier transform, which is trivial because of the delta 
function in the integrand. This yields 
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It is now useful to write explicit expressions for the spin-weighted spherical harmonics 
fl"3| . First, if Z + m is even, then 



K 0) - ( 2 l±l\ 12 [d-m)l(l + m)l] 1/2 

/ml2 ' J_l J V 47r / (Z-m)!!(Z + m)!! ' 



(8) 



while YJ m (§, 0) = if Z + m is odd. Second, if Z + m is odd, then 
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although „iY; m (|, 0) 7^ for I + m even, its expression will not be needed. Finally, we note 
that xY lm (\,G) = (-l) l+m ^Y lm ^,0). 

Looking at this last result and the fact that 1 J = — (_i J), we see that the quantity within 
the square brackets of equation (^) may be simplified as follows: it vanishes identically for 
l+m even, and it is equal to 2_i J _i^ m (f , 0) when l+m is odd. Using this while substituting 
equations (0) and (0) into equation flU), we obtain 
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is the quadrupole-formula result. We also have 
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if I + m is odd. The symmetry r\i_ m = r]i m , used in equation (|I~Ci|), can easily be established 
from the previous expressions. 

Equations (|iTi|)-(|i~4"|) constitute the dressed multipole truncation of the gravitational- wave 
luminosity. The right-hand side of equation ([K]) is an infinite series in powers of v, and our 
next task is to examine the convergence of this series. 

III. CONVERGENCE OF DRESSED MULTIPOLE TRUNCATION 

The velocity dependence of all the terms in equation (|i0|) is made clear by substituting 
equations (0) and (^|), which gives 



e=e q y: 



2(1-2) 



E 



(aim + bimV 2 + c lm v 3 ) I + m even 



1=2 



\ \ V 2 (di m + eimV 2 + flmV 3 ) I + TU odd 



(15) 



The various coefficients are 
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We now express equation (|I5|) in the following form, which displays a more transparent 
grouping of terms: 
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The dressed multipole truncation therefore splits E into three separate series: Si(v), the 
bare multipole series, v 2 S 2 (v), the 1PN correction, and v 3 Ss(v), the 1.5PN correction. It is 
important to point out that in equations (p3|)-(p5|), the sums over m are restricted to values 
such that I + m is even. 

The convergence of the dressed multipole series may be established by analyzing the 
convergence of each of the partial series. We begin with Si(v), the bare multipole series of 
reference ||12|| . Its convergence is established in the following manner. First, it is shown (nu- 
merically) that the terms in the square brackets of equation (^) approach the approximate 
value 1.01 as I — > oo. This implies that the behaviour of the series for large I is dominated 
by the factor in front of the square brackets. The series may therefore be approximated by 
S±(v) ~ v 2 ( l ~ 2 ^aa, which is shown (using the Cauchy ratio test) to converge for 

v < 2/e ~ 0.7358. 



(26) 



The convergence of S 2 (v) is established in the same manner. First, one shows numerically 
that the terms in the square brackets of equation (E3) also approach 1.01 as / — > oo. The 
series is then approximated by S 2 (v) ~ 

and the convergence of this series 
is again analyzed via the Cauchy ratio test. Using equation fllTD , we find that &z+i,z+i/&22 
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is equal to ai+ij+i/au multiplied by a ratio of terms quadratic in I. In the limit / — > oo, 
h+i,i+i/bu ~ a>i+i,i+i/ a u an d we immediately conclude that the convergence of S^v) is also 
determined by equation (p6|). 

The same is true also for S^(v). Again, it can be shown that the terms in the square 
brackets of equation (|25| ) approach 1.01 as I — > oo. The series is then approximated by 
S3(v) ~ v 2 ( l ~ 2 'cu. From equation (|18|) we see that in the limit I — > oo, the ratio of 
successive coefficients is equal to ai+ij+x/au. This allows us to conclude that the convergence 
of S 3 (v) is also determined by equation (|26|). 

In summary, when 1PN and 1.5PN corrections are added to the bare multipole series, 
E splits naturally into three separate series, which all converge for v < 2/e. Therefore, the 
dressed multipole truncation of E must also converge in this interval. We note that this 
convergence test ignores the fact that the combined series S(v ) = Si(v ) +v 2 S2{v) +v 3 Ss(v ) is 
alternating at large I. However, establishing the convergence of the alternating series using 
the Leibnitz criterion does not result in an improvement on equation (^). Furthermore, a 
numerical analysis of the convergence of S(v ) confirms equation (|26|). 



IV. CONCLUDING REMARKS 

We have shown that for circular binaries with small mass ratios, the dressed multipole 
truncation leads to a convergent series for the gravitational- wave luminosity, for all physically 
relevant values of the orbital velocity. It appears likely that this conclusion will remain valid 
for binaries of comparable masses, but a separate calculation is required to confirm this. Our 
results imply that the poor convergence of the true post-Newtonian series does not originate 
from summing over the Newtonian, 1PN, and 1.5PN part of all the multipole moments. 

We see in figure 1 that the dressed-multipole series leads to a significant improvement 
in accuracy, compared with the bare-multipole series. By including the negative 1PN cor- 
rections and the (smaller) positive 1.5PN corrections, our approximation to the luminosity 



coincides (within 1 %) with the numerically "exact" results [10] over the interval v < 0.2. 
One would therefore expect that the incorporation of further corrections would lead to a 
series that converges everywhere in the region of physical interest, and stays accurate over 
a larger interval. 

Meanwhile, the issue of the convergence of the true post-Newtonian series remains an 
open one. 
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FIG. 1. Plots of E/Eq as a function of orbital velocity v. The thick curve is the "exact" 
curve, obtained by numerically integrating the equations of black-hole perturbation theory (see 
reference |l(J). The curve immediately below the exact curve is the 5.5PN curve, obtained from 
the post-Newtonian expansion (carried up to order v 11 beyond the quadrupole formula) of Tanaka 
et al. ||. The first curve above the exact curve is the dressed multipole truncation, obtained 
from equation (|22|). Finally, the highest curve is the bare multipole truncation, obtained also from 
equation (|22|), but with S2(v) and Ss(v) set to zero. We see that the accuracy of the dressed 
multipole truncation is far superior to that of the bare multipole truncation, but that it does not 
match the accuracy of the 5.5PN approximation. 
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